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1. Introduction 



Clusters in intense laser fields were found to absorb laser energy particularly well. The 
ions inside such clusters typically assume higher charge states than expected from a 
single, isolated atom in the same laser pulse. The reason lies in the complex interplay 
between the electric field inside the cluster and ionisation, both "inner ionisation," 
that is, the removal of electrons with respect to individual ions, and "outer ionisation" 
(i.e., removal of electrons from the cluster as a whole). In experiments at wavelengths 

> 248 nm highly energetic electrons H, ions ji El IHl 13 El , photons P UHl HI 1121 , 
and neutrons originating from nuclear fusion ^31 were observed (see Refs. [Ti| ITH] for 
reviews). The first Xe cluster experiments in the VUV regime at high laser intensities 

> lO^^Wcm"^ were performed at the DESY free-electron laser (FEL), Hamburg, 
Germany ^Hl, showing higher charge states than expected from single atoms as well, 
despite the tiny quiver energies of electrons at such short wavelength. 

The interaction of moderately sized clusters with strong laser fields is commonly 
simulated within the framework of molecular dynamics where the inner ionisation is 
incorporated via rate equations or a sequential, single active electron-approach (see, 
e.g., imilHllinillllEIlESEniEll). instead, in this work a semi-classical approach with 
an explicit treatment of the multi-electron inner-atomic dynamics is pursued. In that 
way field or collisionally induced ionisation, recombination, non-sequential ionisation, 
and other inelastic processes are included self-consistently. 

The paper is organised as follows. In section El the semi-classical cluster model is 
introduced before in section Elthe numerical results, especially the ionisation mechanism 
and the energy absorption in clusters at 800 and 100 nm are discussed. We conclude in 
section HI 



2. Semi-classical model 



The rare gas cluster consisting of Na ions of mass M at the positions Ri, 1 < i < Na, 
and Z "active" electrons per ion at the positions Vj, I < j < ZNa in a laser field E(t) 
is modelled by the Hamiltonian 
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(atomic units are used until noted otherwise). Since "classical atoms" with more than 
one electron are generally unstable, a momentum-dependent potential 
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is introduced j2Sl that approximately enforces the Heisenberg uncertainty relation 
when applied in the form VH{fij,pij) = V{fij,Pij,^H,aH,fJ^ei), where Vij = jr^ - Rj\, 
pij = \{Mpi — Pj)/{1 + M)\ are relative position and momentum, respectively, and /lei 
is the reduced mass of electron and ion. Spin effects are neglected but it is worth 
noticing that Pauli-blocking could be incorporated in the model [SEl I2Z1 I2H]- The 
bigger the "hardness parameter" an is chosen, the more severely the uncertainty relation 
"f^ijpij > is fulfilled. However, big an render the equations of motion stiff so that 
in practical calculations a compromise has to be made. C,h is used as a free parameter 
for adjusting the relevant physical properties (e.g., ionisation potentials) of the atomic 
species under study. 

The results presented in this work are obtained by solving the equations of motion 
for all electrons and ions corresponding to the Hamiltonian for both a model Xe 
atom and a model Xe54 cluster. Only the Z = 6 electrons in the 5p valence shell of 
each Xe atom are considered "active." This is a reasonable assumption for the laser 
intensity regime discussed in the present work and significantly reduces the numerical 
effort, an = 2 and C,h = 2.33 was chosen which leads to the six ionisation potentials 
Si = 0.16, 0.55, 1.11, 1.63, 2.22, 2.65 for an isolated Xe atom in this model while the 
real ionisation potentials are (see, e.g., [211 EOl) ^7 = 0.45, 0.77, 1.18, 1.69, 2.09, 2.64. 
The first two ionisation potentials of the model are in poor agreement with the correct 
values. Unfortunately, we found it not possible to tune all six ionisation potentials to the 
correct values by varying the only free parameter (the parameter an hardly affects 
these values j2Zj). Since it is the generation of the higher charge states in clusters that 
is of particular interest to us, C,h was chosen to optimise the higher ionisation potentials. 

3. Results and discussion 
3.1. Single-atom results 

The stationary {r = p = 0) electronic configuration of the single model atom is face- 
centred cubic, each of the six electrons having the same distance \r\ = 1.56 to the ion 
located in the centre of the cube. The canonical momentum is finite, \p\ = 1.49. In 
order to investigate the enhanced (or reduced) ionisation in clusters one needs to know 
the probabilities for the various charge states of the single, isolated atom. In Fig. Q the 
probabilities to find a certain charge state are presented for the two laser wavelengths 
800 and 100 nm. The laser pulse was of trapezoidal shape with an pa 8 fs up and down- 
ramp (corresponding to 3 and 24 laser cycles at 800 and 100 nm, respectively) and a 
T'const = 26.5 fs (10 and 80 cycles, respectively) fiat top. The charge state distributions 
in Fig. n were obtained by simulating an ensemble of 100 randomly rotated single atoms 
in the laser field. Note that the electron dynamics in the laser field and thus the final 
charge state in general depends on the orientation of the initial, classical electronic 
configuration with respect to the laser polarisation Bx- 

The so-called "saturation intensity" for a certain charge state is defined as the laser 
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intensity where the charge state reaches its maximum probabihty before it decreases 
in favour of the next higher charge state. The saturation intensities for charge states 
Z higher than 2+ in Fig.^ agree well with what is expected from the Bethe rule jHI] 
Jsat = Sf/llGZ"^) while the first two charge states appear at higher intensities. This is 
likely due to the incomplete screening of the full nuclear charge by the other electrons 
and a non-negligible down-shift of the binding energy (classical Stark effect) for the 
first two electrons in our model. Note that in the derivation of the Bethe rule a Stark 
shift is not taken into account. However, in real experiments on Xe in 800 nm laser 
light jnij also the saturation intensities for the lowest charge states obey the Bethe rule. 
This means that, by fortunate coincidence, the too low binding energies of the first two 
electrons are partly compensated by incomplete screening and classical Stark effect so 
that the first two saturation intensities in our model are less underestimated as one 
could expect from the ionisation potentials alone. At 100 nm (Fig.^), the maximum 
probabilities to find successive charge states in our model Xe atom lie closer together 
so that broader charge state distributions arise than at 800 nm. 

In Fig. 121 focus-averaged charge state distributions at 800 nm are compared with the 
experimental data in for 200 fs laser pulses. The overall agreement is surprisingly 
good. Note, however, that the experimental laser intensity had to be multiplied by a 
factor 0.4 for matching our model results so that only good qualitative agreement can 
be claimed. 

At 100 nm the energy of a photon is close to the first ionisation potential of real 
Xe. It would be not surprising if our semi-classical model breaks down at such short 
wavelengths where ionisation becomes a single-photon process. Although the model 
possesses eigenfrequencies, these have in general nothing to do with quantum mechanical 
excitation or ionisation energies. First results for 98 nm were obtained at the DESY- 
FEL in Hamburg, Germany (note that in our semi-classical approach it does not 
matter whether one uses 98 or 100 nm wavelength). A comparison suggests that our 
model overestimates the saturation intensities. However, the uncertainty in the peak 
laser intensity and the spiky structure of the FEL pulse inhibits a straightforward 
comparison between theory and experiment. Moreover, it should be emphasised that 
not the absolute values for threshold or saturation laser intensities are the main focus 
of our studies but the relative differences between single atom and cluster ionisation. 

3.2. Cluster results 

Clusters consisting of Na = 54 model Xe atoms in the laser pulse of shape and duration 
as described in the previous subsection were simulated. The initial cluster configuration 
was assembled by attaching randomly rotated single-atom electron configurations to 
the ions. The ion positions Hi were taken from normalised Lennard- Jones cluster 
calculations with each radial vector multiplied by 7.127 so that a Wigner-Seitz 

radius rws ~ 4 was obtained. More details on the actual implementation of our cluster 
model can be found in Ref. [5^ . 
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Figure 1. Probabilities to obtain charge states 1+, 2+, 6+ vs laser peak intensity 
for the Xe model atom in a trapezoidally shaped laser pulse of T « 42 fs duration and 
wavelength 800 nm (a) and 100 nm (b). Multiply the laser intensity by 3.51 x 10^^ to 
obtain the corresponding value in Wcm^^. 




Figure 2. Numerically determined and focus averaged ion yields vs laser intensity (+). 
Dots are experimental data from 32 . The experimental intensities were multiplied by 
a factor of 0.4 . 



3.2.1. lonisation mechanism In order to get a deeper insight into the cluster ionisation 
dynamics it has proven useful to distinguish between three groups of electrons. First, 
there are those electrons that left their parent atom. We consider an electron to 
contribute to this "parent atom ionisation" if it is more than a Wigner-Seitz radius 
away from its parent atom. Electrons in the second group are more than a Wigner-Seitz 
radius away from any ion and thus define the averaged ion charge state 

1 Na ZNa 

= (^iVa - E E ©Ks - \rj - m) (3) 
where 6 is the step-function. 
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Figure 3. Cluster ionisation vs time for tiie peak laser fieldstrengths E = 0.08 and 
A ~ 800 (a) and 100 nm (b). The average number of electrons per ion that left their 
parent atom ("parent"), any atom ("charge state"), or the entire cluster ("outer") are 
presented. Vertical dashed lines indicate the end of the laser pulse. Horizontal green- 
dotted lines mark the charge state of an isolated atom after the same laser pulse. See 
text for discussion. 



Finally, there are those electrons that are far away from the entire cluster. We 
count these electrons by calculating 

ZNa 

iVouter = Yl ~ max{i?i} - rws) (4) 
i=i 

which equals the total charge of the cluster that is left behind. The charged cluster 
is unstable and Coulomb-explodes. Under the conditions investigated in this work, a 
typical time within which the cluster doubles its radius is several tens of femtoseconds. 
In the definition (jH) the cluster radius is taken as -Rciuster = max{i?j} + r-yvs where 
max{i?j} is the distance from the outermost ion to the cluster centre. 

In Fig. inithe parent atom ionisation, the average charge state, and outer ionisation 
are plotted vs time for the laser intensity I = E'^ = 6.4 x 10~^, corresponding to 
2.2 X lO^^Wcm"^, and the two wavelengths A = 800 and 100 nm. Common features in 
both results are the fast removal of all active electrons from their parent atom and the 
moderate outer ionisation 2 electrons per ion are removed from the entire cluster). 
The removal of all six electrons from their parent atom might be surprising in view 
of the fact that single, isolated atoms in the same field loose not even three electrons 
(see horizontal, green-dotted line in Fig. IHK)- At lower laser intensities this discrepancy 
between single atom ionisation and parent atom ionisation in clusters is even more 
pronounced. In fact, we found that it is in general sufficient that the single, isolated 
atom looses one electron to trigger complete parent atom ionisation in the cluster. 
However, parent atom ionisation is not an observable! Electrons are indistinguishable, 
and, indeed, other electrons are trapped by the ion. The average charge states — 
depending on the balance between ionisation and recombination — are significantly lower 
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in Fig. 121 Looking at the particle trajectories one infers that the electrons hop from one 
ion to the other, undergoing inelastic collisions and getting temporarily trapped. This 
is expected because a low temperature, strongly coupled plasma of electron density Ue 
and Debye-length is formed whose plasma parameter l/{ne\%) is close to unity or 
even higher. Consequently, a Debye-sphere is populated by only a few electrons, and 
the potential energy of the electrons is comparable to their kinetic energy. 

One may object that complete parent atom ionisation suggests that more active 
electrons per atom should be incorporated in the simulation in the first place. However, 
what is important is that the charge states 



remain smaller than the highest allowed charge state Z. Whether an electron screens 
(or is bound to) its parent atom or another one is insignificant. Hence, increasing the 
number of active electrons per atom while < Z introduces an unnecessary large 
overhead of bound electrons to be propagated. The stronger coupled the plasma is, the 
greater the difference between the "parent atom ionisation degree" and the charge states 
(0) is expected. 

The average charge state Zavcs (indicated by "charge state" in Fig.E)) was calculated 
geometrically by counting all electrons farther away than r-ws = 4 from any ion. This 
geometric concept is meaningful only when the cluster has sufficiently expanded so that 
it becomes unlikely that free electrons occupy accidentally the test spheres of radius 
r-ws around the ions. In fact, the cluster expansion is the reason why the charge states 
appear to increase even after the end of the laser pulse (indicated by the dashed vertical 
lines). Hence, only the asymptotic values of the "charge state"-curves in Fig. El are the 
expected average charge states that would be measured in a real experiment. In other 
words: when the ions are well separated (as it is the case in an experiment, shortly before 
they hit the detector) all electrons within a sphere of radius r-ws around the ions can be 
safely considered bound so that © is indeed the charge state that would be measured, 
and, consequently, ^avcs = J2 

Z^J^/Na would be the average charge state provided 
that all ion species are detected with the same efficiency. In order to figure out whether 
the cluster ionises more efficiently than the single, isolated atom, the average charge 
state Zavcs can be compared with the corresponding results presented in the previous 
subsection. These single atom results are included in Fig. El (horizontal, green-dotted 
lines) . 

It is seen in Fig. IHl that, for the particular laser intensity chosen, at 100 nm 
the cluster ionisation is less efficient than the single-atom ionisation. The "relative 
ionisation," that is, the average charge state in the cluster divided by the corresponding 
single atom charge state, is shown in Fig. |31 as a function of the laser intensity. A value 
above unity indicates that ionisation in the cluster is more efficient than ionisation of 
the isolated atom. At both laser wavelengths the relative ionisation in the cluster is 
particularly high when the first free electrons are generated. The ionisation rate for the 
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Figure 4. Cluster ionisation normalised to the single, isolated atom ionisation (see 
inset) vs laser intensity for the two wavelengths A = 800 (red) and 100 nm (blue). 
Average charge states (solid lines) and outer ionisation per cluster atom (dashed) are 
shown. See text for discussion. 



very early electrons in the cluster necessarily equals the one for isolated atoms since 
there is not yet the potential from the other ions. However, as soon as these very early 
electrons leave, the cluster charges up, leading to strong electric fields, especially at 
the cluster boundary, that enhance further inner ionisation ("ionisation ignition" jl7j . 
details follow below). This is the reason for the strong relative ionisation at intensities 
/ < 6 X 10~^ for both wavelengths. 

At higher laser intensities inner ionisation increases but the electrons are less 
efficiently transported away from the cluster as a whole — especially at short wavelengths 
because the quiver amplitude of the freed electrons is low then so that outer ionisation 
can proceed through thermionic emission only. The electron temperature that is reached 
is sensitive to the pulse duration. For the 42fs-pulses and the laser intensities in the 
present work the electron temperatures at 100 nm always remain below 1.7 atomic units. 
During the cluster expansion, the electrons trapped by the cluster potential are cooled, 
and recombination leads to only moderate final charge states (the higher, transient 
charge states should be experimentally accessible by measuring the emitted X-rays). As 
a result, the relative ionisation at 100 nm drops below unity, that is, the atoms inside the 
cluster ionise — on average — less efficiently than an isolated atom does. At 800 nm and 
the same laser intensity relative ionisation remains above unity, mainly because single 
atom ionisation is lower (see inset in Fig. Ej) and the quiver energy suffices to free the 
electrons from the cluster potential (see red, dotted outer ionisation- curve in Fig. H}. 
Finally, at sufficiently high laser intensities all active electrons are removed from both 
the cluster atoms and the isolated atom so that relative ionisation approaches unity. 
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In order to analyse the origin of the efficient ionisation inside the cluster in more 
detail it is useful to look at the external forces a bound electron experiences. This is 
done in Fig. El for the two cases of an electron bound to an ion at the cluster boundary 
and to an ion sitting in the cluster centre. The laser parameters were as in Fig. Eh- 
The forces experienced by the electron are divided into a macro field-part which is 
made up by the laser and all electrons and ions that are outside a test sphere of radius 
'"ws = 4 around the ion of interest, and into a micro field-part originating from electrons 
passing through the test sphere within tth = 8 atomic time units. The threshold time 
tth is introduced in order to estimate the importance of coUisional ionisation of the type 
e -|- y4^+ — > 2e + ^4^^+^^+ without allowing processes where the incoming electron is 
captured to affect the micro field (like e -|- A^~^ — > e' + A^~^, for instance). Other ions 
never pass through the test sphere because the ions do not overtake each other (see 
below for an explanation). 

From Fig. |3] one infers that for the laser parameters used in Fig. / = = 
6.4x 10~^ and A = 800 nm, the average cluster ionisation is more efficient than in isolated 
atoms. Fig. EJd explains why. During the ffist four laser cycles the test atom exclusively 
"sees" the laser field because ionisation is still very low so that there are no free electrons 
around. Then the ffist electrons are removed (see the parent atom ionisation in Fig. Efc) 
and the macro field significantly exceeds the value one would expect from the laser field 
alone. Looking at the centre of mass in laser field direction (Fig. Ek) of all electrons 
inside the cluster, it is seen that the pronounced maxima in the macro field occur at 
times where most electrons are on the opposite side of the cluster so that the screening 
of the positive cluster charge is low (see dashed lines to Fig.|3^). Half a laser cycle later, 
when most electrons move to the side of the cluster where the test atom is located, the 
macro field drops dramatically. This clearly shows that a dynamical version of ionisation 
ignition should be adopted, as proposed in Ref. jSZ]. As outer ionisation increases, the 
maxima of the macro field increase as well. This is due to the positive ion background 
that pulls electrons from the cluster boundary into the interior, i.e., the "standard" 
ionisation ignition mechanism. The force -Fignition = A^outer/-R^, where R is the radial 
position of the test atom under consideration and iVouter is the cluster charge as defined 
in (^, is included in plot (b). It is seen that the overall temporal evolution of the macro 
field is in good agreement with Fignition- Deviations, especially at late times, are due 
to the inappropriate assumption of a homogeneously charged sphere for an expanding 
cluster consisting of only 54 ions. The micro field, caused by electrons passing through 
the test sphere, does not contribute to a net increase of the charge state but only to 
fluctuations (see Fig. Eh)- 

The test atom in the centre clearly experiences no ionisation ignition. Although 
this atom looses all its "initial" electrons (see "parent" -curve in Fig. EF) the overall 
charge state remains low, in fact, lower than for an isolated atom in the same laser field. 
This shows again that calculating charge states from parent atom ionisation would lead 
to wrong results. The micro field is stronger in the cluster centre because the electron 
density is higher there. However, the micro field fluctuations again do not increase the 
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Figure 5. Panels (b) and (e) illustrate the absolute values of the external forces 
vs time, as seen by an electron bound to an ion at the cluster boundary and in the 
cluster centre, respectively. The electric field envelope of the laser pulse {E — 0.08, 
A — 800 nm, marked "pulse envelope"), the force exerted by all electrons and ions 
outside the test sphere of radius r-yvg = 4 plus the force of the laser ("macro"), and the 
force exerted by all electrons passing the test sphere within a time ith = 8 ("micro") 
are shown. In (b) the "ignition" force expected due to the positively charged cluster is 
included as well. The centre of mass of all electrons inside the cluster Xcm is plotted 
in (a) and (d). Parent atom ionisation and charge state of the cluster atom under 
consideration are shown in panels (c) and (f). 



net charge state on longer time scales. Hence, coUisional ionisation is not responsible for 
the increased ionisation at 800 nm in Fig. H] neither for atoms at the cluster boundary 
nor for atoms in the centre. 

From Fig.|3]one infers that at the shorter wave length A = 100 nm and / = 6.4x 10~^ 
the average charge state inside the cluster is slightly less than for the isolated atom. In 
Fig- El the external forces experienced by a bound electron at the cluster boundary and 
in the centre are presented for this case. Because of the higher laser frequency the 
centre of mass-excursion is less pronounced so that collective electron dynamics inside 
the cluster does not build up. Note that the centre of mass-excursion is actually less 
when the laser is on than it is after the laser pulse, as it is expected for a driven harmonic 
oscillator when the driving frequency is higher than the eigenfrequency (i.e., the Mie 
plasma frequency in our case). Outer ionisation during the laser pulse proceeds through 
collisional heating (see section 13. 2. 2l below) so that the cluster charges up and the macro 
field again nicely follows -Fignition in Fig. Eb- The final charge state is maximal, i.e., 6, 
but one has to keep in mind that also the isolated atom looses already ~ 5 electrons in 
this laser pulse. The number of electrons inside the test sphere around the ion in the 
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Figure 6. Same as in Fig. but for A = 100 nm. Note tliat tlie time interval siiown 
is twice as long now. 



cluster centre strongly fluctuates, slightly drops (corresponding to an increase of Z^g in 
Fig. IHT) after the laser pulse because of the previously mentioned localisation effect of 
the laser and cluster expansion, but remains below the isolated atom result. 

The fact that the charge states vary from high (at the cluster boundary) to low 
(in the cluster centre) already suggests that the charge state distribution should be 
much broader than in the single, isolated atom-case. This is conflrmed in Fig. [7| where 
the charge state distributions for both laser wavelengths can be compared with each 
other and the corresponding single, isolated atom-results. At the shorter wavelength 
the average charge state in the cluster falls behind the isolated atom ionisation because 
the electron temperatures reached are not high enough to yield sufficient outer ionisation 
through thermionic emission. As a result, the electron density inside the cluster 
potential remains high and, through recombination, leads to a flnally broad charge 
state distribution with a mean value lower than in the single atom-case. 

Another consequence of the fact that the charge states increase from the cluster 
centre towards the cluster boundary is that during the Coulomb explosion the ions do 
not overtake. In the deflnition of the micro fleld above ions therefore needed not be 
taken into account. 

From Fig. E]it can be inferred that in our model only for / < 0.002 the average 
cluster ionisation is more efficient than the isolated atom ionisation. It might be this 
regime into which the enhanced ionisation measured in Xe clusters at 98 nm for laser 
intensities up to 7 x 10^'^ Wcm~^ falls into. It would be interesting to verify 
experimentally whether at higher laser intensities relative ionisation in small clusters 
really decreases, as it is suggested by our results in Fig. HI for J 0.006. 
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Figure 7. Charge state distributions a,t I = E'^ = 6.4 x 10~^, A = 800 (red) and 
100 nm (blue) for the single, isolated atom (a) and the cluster consisting of 54 atoms 
(b) . The mean values are indicated by the narrow bars at the top of each panel. 



3.2.2. Energy absorption The energy absorbed by the cluster in the laser field E{t) is 

Na ZNa 

AS = n{T) - n{o) = J2As,+ J2 (6) 

i=i j=i 

where T is the pulse duration, and 

AS, = Z j^E{t) ■ iti dt, ASj = - j^E{t) ■ Vj dt (7) 

are the absorbed energies of the individual ions and electrons, respectively. Hence it is 
necessary for laser energy to be absorbed that the total current J = Z J2i Ri ^ J2j 'f'j 
acquires a phase lag 7^ ±n/2 with respect to the driving field E(t) so that AS = 
Jq E{t)-Jdt ^ 0. There are several mechanisms that can lead to dephasing: electron-ion 
collisions (i.e., inverse bremsstrahlung) , electrons "colliding" with the cluster boundary 
^JEH]! and electrons colliding with the cluster as a whole [33 EDI- In the following the 
relative importance of these mechanisms is analysed for our case of a small cluster in a 
low- and high-frequency laser field. 

Fig. IHl shows the total absorbed energy © vs the peak field strength of the laser. 
Absorption is more efficient for the short wavelength up to -E ~ 0.13. This is because 
inner ionisation is stronger for the high frequency laser field so that both electron density 
and average charge states are higher than at A = 800 nm. However, the long wavelength- 
result overtakes around E ^ 0.13. This can be attributed to the increasing importance 
of electron-cluster and electron-cluster boundary collisions, as it is evidenced in the inset 
in Fig. ISl where the ratio ASoa/AS with 

ZNa Na 



AS. 



I m ■ f, li e(|r, - R,\ - rws) dt (8) 
•^0 7=1 j=i 
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Figure 8. Laser energy absorbed by the Xe54 model cluster for the two laser 
wavelengths vs the laser field amplitude E. The fraction A£oa,/A£, contributed by 
electrons while they are farther away than r^g = 4 from any ion ("outer absorption"), 
is shown in the inset. 



is plotted. The "outer absorption" ASoa measures the laser energy that is absorbed 
by mechanisms different from standard inverse bremsstrahlung through electron-ion 
collisions because only electrons while they are farther away than rws = 4 from any 
ion contribute in (jH)). Since the maximum excursion of free electrons in a laser field 
is ~ E/uj"^ it is expected that cluster boundary-effects are more important at 800 nm. 
In fact, even at the highest field amplitude shown in the inset of Fig. IHlthe excursion 
EjixP' = 1.4 at 100 nm is still much smaller than the initial cluster radius 15.4 so that it 
is clearly the thermal velocity that dominates the dynamics at 100 nm. 

Also connected to the more efficient absorption at 800 nm in stronger fields is the 
fact that freed electrons are rescattered from the entire cluster while they oscillate with 
large amplitude in the laser field. At 100 nm instead the electrons leave the cluster 
through thermionic emission, and the quiver amplitude in the laser field is too small for 
driving them back to the cluster. Consequently, the freed electrons drift away and are 
thus lost for further energy absorption. 

Fig. ini illustrates with sample trajectories when and where electrons absorb 
predominantly laser energy. The absorption rate averaged over a laser cycle T, 

^,(t) = -- / E{t').r,dt\ (9) 

1 Jt-T/2 

is a measure for the efficiency of laser energy absorption and is used to colour-code the 
trajectories and curves in Fig. 121 

In the long wavelength case collisions with the cluster as a whole clearly dominate. 
The biggest jump in the absorbed energy in Fig. IHb occurs between points A and B 
where the electron crosses rapidly the cluster and large-angle electron-ion collisions are 
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Figure 9. Trajectory (a,c) and absorbed energy vs time (b,d) of the cluster electron 
that contributes most to the total energy absorbed by the cluster at 800 (a,b) and 
100 nm (c,d) for a laser field amplitude E — 0.07. Colour indicates whether the cycle- 
averaged absorption rate is low (dark colours) or high (light colours). The "outer 
absorption" contribution (dashed) is also included in plots (b) and (d). For the sake 
of clarity, the electron trajectories in (a) and (c) are shown only for times i G [0, 1450] 
and t e [400, 1000], respectively. 



obviously absent. Consequently, the fraction of "outer absorption" is particularly high 
(see dashed curve in Fig. Eb)- Note, that the time the electron takes to pass the cluster 
from A to B is of the order of a laser period. Hence, the so-called impact approximation 
(collision time -C laser period) does not apply, and efficient energy absorption is possible 
also in small angle scattering events. 

The electron trajectory in the short wavelength case (Fig. Efc) is rather erratic due 
to many individual electron-ion collisions. The electron temperature Tg towards the 
end of the laser pulse is ~ 1.3 (corresponding to ~ 35 eV) while the quiver energy in 
the laser field t/p is only 5.9 x 10"^. Consequently, the motion of the electron inside 
the cluster is dominated by the thermal velocity rather than by the laser field. For the 
long wavelength instead, the free electron excursion is already greater than the (initial) 
cluster radius, the quiver energy Up is 0.38, and the final electron temperature is of the 
same order of magnitude, namely Tg ~ 0.7. 

3.2.3. Collision frequency For making comparison with analytical theories of inverse 
bremsstrahlung it is now attempted to calculate a time-averaged collision frequency 
(z/) from the numerically determined total cluster energy £{t). Absorption rate £ and 
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2000 



Time (a.u.) Time (a.u.) 



Figure 10. Total cluster energy vs time for E — 0.3 at 800 (a) and 100 nm (b). 



collision frequency are related through (see, e.g., jUj) 

^^ = {r-E) = {ur')^2{u)U,, (10) 

where (■ ■ ■) indicates time-averaging over one laser cycle, is the number of electrons 
undergoing collisions, and Up = E'^/{Auo'^) is the quiver (or ponderomotive) energy. 
Calculating the collision frequency iy) is complicated by the fact that both £ and N,, 
vary in time. The total cluster energy as a function of time at 800 and 100 nm is 
shown in Fig. El for the case E = 0.3 . It is obviously not possible to determine a 
constant absorption rate S from the overall slope of these curves which would be valid 
throughout the laser pulse. In particular at 100 nm the energy absorption strongly 
saturates because almost all active electrons leave the cluster. At low laser intensities, 
shortly after ionisation sets in at all, the behaviour is the opposite; then S increases 
during the pulse. 

Because of these difficulties related to the time-dependent electron density inside 
the cluster, an averaged single-electron absorption rate was introduced as 

AS , , 

where Tconst is the duration of the constant-intensity part of the trapezoidal laser pulse. 
For the number of electrons, N^, = NaZ^^^siT) with T the total pulse duration, was 
chosen at 800 nm and iVg = NaZ^^csiT) — NonteriT) at 100 nm because free electrons may 
be driven back to the cluster in the former case while it is less likely in the latter. Taking 
the final values for Z^vcsit) and Nonterit) may appear a too rough estimate. However, 
the average charge states, ion densities, and electron temperatures to be inserted in the 
analytical collision frequency formulas were also taken at time T so that consistency is 
maintained. 

We compare the numerically determined collision frequency 
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with 

3v^ 



\Ki ) = TT-^ In A (13) 



where rii is the ion density, and f is a typical electron velocity, v = (3Te + /u'^Y^'^, 
which might be dominated by either the electron temperature Tg or the driving laser 
field of amplitude E. For the Coulomb logarithm either 02| 



or 

V + v' 



In A/ = In 



V — v 



(fast e) (15) 



with C the Euler constant and v' = (v'^ + 2ujY^'^ was taken. The expression for fast 
electrons is expected to hold for v ^ ^avcs and u <^ v"^, while the one for slow electrons 
should be valid for v <^ 2'avcs) ^ ^^/^avcs 

Since the condition a; <^ f ^ is in general not fulfilled at 100 nm, we also make 
comparison with the expression derived in [32] I 

, (K). ^ 87r3/^n. (2Zl\'^' £(1/3) 

\ / ^5 35/6(2re)l/2 / r(2/3) ^ ' 

which is supposed to be valid as long as ^ t>^/Zavcs- 

The result is shown in Fig. ^2 It is seen that at 800 nm none of the electron- 
ion collision frequencies (fT!?|) - (fTT)|) agrees with the numerical result, especially around 
E ^ 0.06 where both electron temperature and quiver energy are relatively low (0.6 
and 0.28, respectively). If the collision frequency (z/gj) was as high as predicted by these 
formulas, electron-ion collisions should clearly exceed the energy absorption by electron- 
cluster boundary collisions which, however, is not the discussed in section IH. 2. 21 

The numerical result at 100 nm is in good agreement with (I'l.f) for not too low laser 
intensities whereas (z^ef underestimates the numerically determined collision frequency. 
In Ref. pll formula ()16|1 was found to overestimate coUisional absorption which was 
attributed to screening effects not included in the derivation of ()16|). In fact, the Debye- 
length Xd = [Te/ {AnneY^'^ is of the same size or even smaller than the interionic distance 

1 /3 

n~ so that one would expect (fTBj) to overestimate the collision frequency — which it 
does not in our case. 

The reason for the discrepancies between the simulation results for the collision 
frequency and the analytical values is not clear (at least to the author). Only if one 
allows for an unreasonably generous adjustment of the effective charge states (commonly 
justified with screening of the nuclear charge) acceptable agreement could be achieved in 
the 800 nm case. Also unclear is why, at 100 nm, eqs. ()13pi4j) work better than expected 
while the seemingly more appropriate formula (fT^ gives poor agreement. 

X Note that Eq. (49) in ,42, differs from our Eq. H10|l by a factor of 2. This translates to the expression 
for (^ff V 
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Figure 11. Collision frequencies vs peak field strength of the laser at 800 (a) and 



100 nm (b). Numerically determined collision frequency (solid, bold), 
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In Ref. ^42, enhanced laser energy absorption was attributed to electrons scattering 
on an effective potential that changes from Coulombic with the full nuclear charge in 
the interior to a Debye-screened potential with the actual charge state in the exterior. 
Electrons with small impact parameter thus "see" a nuclear charge higher than the 
actual charge state. This effect, of course, is not specific to laser-cluster interaction but 
is a general issue in plasma physics. It is automatically included — may it be important 
or not — in our numerical calculations (at least as long as the number of actively treated 
electrons is sufficient) while it is not incorporated in molecular dynamics simulations 
dealing with rates for inner ionisation and effective charge states only. 

In our work, ionisation ignition was found to be the dominant ionisation mechanism 
in clusters also at short wavelengths. This agrees with results presented in |24j . 
Ionisation ignition, upon relying on outer ionisation, is a finite-size effect not considered 
in ^3] • It is the mechanism for outer ionisation that depends on the laser wavelength and 
ranges from collective electron dynamics (long wavelength) to collisional absorption and 
thermionic emission (short wavelength). Hence, inverse bremsstrahlung only indirectly, 
via outer ionisation, leads to the generation of high charge states in clusters. 

4. Conclusions 

Results from numerical simulations of a semi-classical Xe54 model cluster in short and 
intense laser pulses at 800 and 100 nm wavelength were presented and compared with 
the corresponding single atom results. The inner-atomic dynamics was treated explicitly 
so that no rates for field ionisation, collisional ionisation, and recombination had to be 
employed but all physical effects, although classical ones only, were included exactly. 

It was found that (dynamical) ionisation ignition is the by far dominating 
mechanism behind enhanced ionisation of atoms inside clusters at both long and short 
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wavelength. The outer ionisation mechanism depends on the laser wavelength. At 
long wavelengths, when the electron excursion is comparable to the cluster radius, laser 
energy absorption proceeds mainly through electron-cluster (boundary) collisions while 
at short wavelengths, where the electron dynamics is dominated by the thermal rather 
than the quiver velocity, inverse bremsstrahlung prevails. 

Collisional ionisation is insignificant. However, under the conditions studied in this 
paper the laser field turns the cluster into a strongly coupled plasma where inelastic 
processes such as electron capture and reemission occur with high probability. As a 
consequence, the charge states remained moderate although all active electrons in the 
simulation left their parent atom. 

The charge state distributions of the cluster ions are broader, ranging from the 
highest charge state generally produced at the cluster boundary down to the lowest one 
produced in the cluster centre. At short wavelengths the small electron quiver amplitude 
and the moderate temperatures reached in short pulses hamper outer ionisation. As 
a result, after cluster expansion the final, mean charge state may lag behind the 
corresponding single atom result. 

Recent work employing the "particle-in-cell" simulation technique [H] indicates 
that at long wavelengths dynamical ionisation ignition (dubbed "polarisation enhanced 
ionisation" there) and electron-cluster boundary collisions (called "laser dephasing 
heating" in ^3]) remain to be the dominant ionisation and absorption mechanisms 
in big clusters of several 10^ atoms as well. 

Finally, the question arises whether the explicit treatment of the inner-atomic 
dynamics leads to results different from those obtained with simpler approaches. On one 
hand, the method used in this work, although semi-classical only, is fully self-consistent 
and clearly superior to approaches where field ionisation, recombination, collisional 
ionisation, etc. are treated through rate equations. Accurate rates for these processes in 
non-equilibrium laser plasmas are hardly available, and the different processes cannot 
even be properly separated. On the other hand, there is a type of simulations where 
the outermost bound electron is considered explicitly while the deeper bound ones are 
accounted for by an effective nuclear charge of the parent atom. The parent atom 
charge state is increased when the electron leaves, and the same procedure starts for 
the next bound electron. These kind of simulations which enforce "purely sequential 
parent atom ionisation" (PSPAI) should lead to similar results as compared to ours 
since non-sequential ionisation does not play a major role. In particular, there should 
be the same important difference between the physically meaningful charge state (0) 
and the non-observable parent atom ionisation. 

If very high charge states are created the method used in the present work becomes 
more rapidly inefficient than those neglecting the multi-electron, inner-atomic dynamics 
because too many electrons per atom need to be simulated from the beginning. In 
such cases the "PSPAI approach" described in the previous paragraph is the method 
of choice, possibly amended to account for explicit recombination so that the effective 
nuclear charges and the number of electrons in the simulation will eventually decrease 
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while the cluster expands. 
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